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Abstract 

We study the transition between Zeeman levels of an arbitrary spin placed into a regular time- 
dependent magnetic field and a random field with the Gaussian distribution. One component of 
the regular field changes its sign at some moment of time, whereas another component does not 
change substantially. The noise is assumed to be fast. In this assumption we find analytically the 
ensemble average of the spin density matrix and its fluctuations. 

PACS numbers: 23.23. +x, 56.65.Dy 



I. INTRODUCTION 



Nanomagnets Fe$ and Mnu attracted much attention last decade. It was shown ex- 
perimentally that they realize a single molecule quantum hysteresis 1 . Some features of the 
hysteresis have been predicted in many theoretical works^i. In particular, it is commonly 
accepted that plateaus on the hysteresis loop are due to the Landau-Zener (LZ) transitions^ 
at avoided crossing of Zeeman-split spin levels in the crystal field of molecular environment. 
This point of view is strongly supported by the experimental observation^ of theoretically 
predicted 8 oscillations of the transition matrix elements vs. magnetic field applied in the 
hard direction. However, it was indicated earlier— that the hyperfine interaction with the 
numerous nuclei in the molecule is comparatively large and can remarkably complicate a 
simple LZ picture of transition. The reason why it did not happen at a specific experimen- 
tal setup^ was recently explained theoretically^. It is indicated in the same work that at a 
smaller rate of the field sweeping the violation of the simple LZ formula seem to be unavoid- 
able. The nuclear spins provide not only a random static field which changes locally Zeeman 
splitting and create a non-orthogonality. Their fluctuations at a little higher temperature 
become fast and can be considered as a random time-dependent magnetic field acting on the 
spin. Another source of the noise is the interaction with phonons 10 , which produces random, 
time-dependent anisotropy. 

Recently several nanomagnets with cubic or almost cubic symmetry were synthesized. 
One of them MoeMng 11 has spin 51/2 (this value can be comparatively easy varied by simple 
substitutions). Another nanomagnet [(trifos) Re (CN) 3 ] 4 [CoCl]^, which we will abbreviate 
as Re4Co4, has a smaller spin 5 = 6. So far no effects of anisotropy were found in magnetic 
measurements even at temperature about IK. It shows that the cubic anisotropy is weak, as 
it could be expected. Thus, the cubic nanomagnets may be the best realization of a large free 
spin. Being placed into a varying magnetic field B z (t) along z— axis and a small constant 
field B x in x— direction, such a spin performs quantum transitions between the states of the 
Zeeman multiplet. The essential difference between this situation and the standard Landau- 
Zener (LZ) problem^ is that at B z (t) = all 2S + 1 states of the Zeeman multiplet cross 
simultaneously, whereas Landau and Zener considered only two-level crossing equivalent to 
S = 1/2. An extension of the LZ theory to the case of higher spins was proposed by Hioe^. 
He has proved that the transition probabilities depend only on S and the dimensionless LZ 
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parameter 

7 = y/gfi B B*/4hB z , (1) 

where g is the Lande factor, is the Bohr's magneton and B z is the time derivative 
(sweeping rate) of the magnetic field taken at the diabatic level crossing, i.e. at the moment 
when B z (t) turns into zero. The transition probabilities display very interesting oscillations 
vs. the parameter 7 and difference of projections m and ml of the orbital moment in the 
initial and final states. These oscillations, which were not analyzed in the work 1 ^, originate 
from the quantum interference of different Feynman paths leading from m to m', which does 
not exist in the genuine LZ problem 14 . Thus, the cubic nanomagnets are simple, but non- 
trivial physical objects. The quantum tunneling in these objects is worthwhile of studying. 
As it was shown for other nanomagnets, the interaction with nuclear spins is essential for 
the dynamics. Though the nuclear relaxation times are sufficiently long in microscopic scale 
(typically milliseconds), the characteristic time of the LZ process r^z ~ B x /B z may be even 
longer. Then the field of nuclear spins can be considered as a fast gaussian noise. In the 
opposite case this field should be considered as a quasistatic random field. Such a case was 
studied theoretically by Prokof 'ev and Stamp 4 -. 

The purpose of our article is to study the influence of the fast noise onto probabilities 
of transitions between the states of Zeeman multiplet for a free spin S in the presence of a 
regular time-dependent field, the same as that in the LZ or Hioe problem. A special case of 
this problem for S = 1/2, when the regular field B x is zero and transitions are completely 
determined by non-diagonal elements of the random field, was solved earlier by Kayanuma 1 ^, 
who has found average values of the transition probabilities. In our previous worki& we have 
solved the same problem for regular and random fields acting together. In this work we 
extend our results for higher spins, find the average non-diagonal matrix elements of the 
density matrix and calculate the fluctuations of all these values, which occur to be strong. 
These calculations became possible due to high symmetry of the problem. A proper group- 
theoretical treatment allows to deal with standard objects, which we call Bloch tensors. 
They are a generalization of the well-known Bloch vector for spin 1/2 problem. 

The plan of the article is as follows. In section 2 we remind the Hioe solution and analyze 
in some details the oscillations of the transition probabilities. In section 3 we formulate 
equations for the density matrix and reduce them to equations for the Bloch tensors, which 
we define in the same section. In section 4 we consider the fast noise acting on a two-level 
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system. Though this problem was considered in our previous article 16 , it is important to 
give a simple analysis of the time scales and to calculate average transition matrix and 
fluctuations, which will serve as a basis for the spin-S* problem. Section 5 is dedicated to 
the solution of the same problem for general spin. In section 6, returning to 2- level system 
(spin 1/2), we discuss the limit of a strong regular field B x , so that the regular transition 
proceeds adiabatically. The noise remains fast in the scale tlz, but its relaxation frequency 
may be much larger than B x . The last 7-th section contains our conclusions. 



II. SPIN S LANDAU-ZENER PROBLEM 



The Hamiltonian of a free spin S with the maximal projection S in an external time- 
dependent magnetic field B (t) reads: 

H s (t) = -Sb(t) (2) 

where b (t) = gfiB^{t)- The key observatio n 13 ! 14 is that this Hamiltonian is an opera- 
tor of time-dependent infinitesimal rotation. Therefore the corresponding evolution matrix 
Us(t,t ) = Texp \—iJ Hs(t')dt'\ is an operator of rotation belonging to the group 50(3) 
acting in its irreducible representation labeled by an integer or half-integer S. Since the 
composition law is the same for any irreducible representation, the resulting evolution op- 
erator represents the same rotation for any spin. The group theory allows to construct this 
matrix for an arbitrary spin if it is known for spin 1/2 (see^, ch. VIII). The multi-spinor 
technique is most appropriate for this purpose. The spin S state can be represented as a 
direct symmetrized product of 2S spin 1/2 states: 



\S, m) 



(S + m)\(S — m)\ . . . 

(|+ + ... + -- ...-) + |+ + ... - +...-) + ...) (3) 



\ (2^)! 

where each ket contains S + m spins up (+) and S — m spins down (— ) and all permutations 
of up and down are performed. Let the SU(2) matrix rotating spin 1/2 states be: 



u 



(4) 



V 



a b 
-b* a* 

with the constraint |a| 2 + \b\ 2 = 1 imposed. The transformation for the state of the spin 
S can be obtained as the direct product of transformations (JH): 



\S,m) -> J < s+ '^ s y - m> ' a s +"'(-b') s - m \S, S) + 

, (s -Si^ s+ '"(-*>*) s -™-v + (-f)^-t) 1) ... (5) 

A general matrix element of the rotation operator (m| Us \m') for the spin S is expressed in 
terms of a, b, a*, b* in the following wa y 17 i 18 : 



(m\ Us \m') 



(S + m)\(S m)\ m 1 +mim' -m nm' -m,m' +m ( r)\ \2 i \ / K \ 

where P^ ,b (x) are the Jacobi polynomials^. The matrix elements possess the follow- 
ing symmetry properties: {—m\Us\—rn!) = (— 1)H+I m 'l (m\ Us |wi')*> K m l % I 771 ')! = 
| (m'| C/s |m)| = |(— m| £/g |— m')\. Equation © displays oscillations of the matrix element 
when the argument 2|a| 2 — 1 varies from —1 to +1. These oscillations are associated with the 
oscillatory behavior of the Jacobi Polynomial. For the number of nodes N(m,m', S) of the 
matrix element (m\ Us \m!) a simple equation is valid: N(m,m', S) = S — max(|m|, \m'\). 
The central matrix element with m = m! = for integer S and \m\ = \m'\ = 1/2 for 
half-integer S has maximal number of nodes equal to S and S — 1/2, respectively. 

Let us specify the problem considering only a narrow vicinity of the diabatic levels crossing 
point, which we accept for t = 0. It is possible if the interval of time tlz, during which the 
transitions presumably proceed is much less than the characteristic time of the field variation 



tn 



B z /B z . The characteristic time of transition can be identified as tlz = B x /B 



Thus, the requirement of short transition time can be rewritten as B x <^ {^z) I 
this requirement is satisfied, one can approximate with high accuracy the magnetic field by 
a linear function of time B z (t) = B z t. 

Landau and Zenei^ have solved such a problem for spin 1/2. In particular Zener have 
determined matrix elements a, b in terms of the Weber function D_ iry 2[e m ^\JCot) for an 
arbitrary moment of time. For simplicity we will focus on the values of these parameters for 
transition from t = —00 to t = 00. According to the Landau-Zener solution 

a = exp(-7r 7 j; b = ^(-if 1 ) (7) 

When 7 varies from to 00, the modulus \a\ changes from 1 to and the argument of the 
Jacobi polynomial in equation © varies from 1 to —1. Using the expression one finds 
the solution for an arbitrary spin S in terms of coefficients a and b, or equivalently, in terms 



B, 



If 



of the Landau-Zener parameter 7. Thus, the transition amplitudes oscillate as function of 
7. The physical reason of these oscillations is the interference between different Feynman 
paths leading from m to m! . For illustration we show here corresponding matrices for spins 
1, 3/2 and 2: 



( a 2 -V2ab -b 2 \ 
V2ab* 2\a\ 2 -l V2a*b 

-b 



<^ 2 -V2a*b* a* 2 



( 



U3/2 



y/3a 2 b 



V3ab 2 



b 3 



\ 



V 



/ 



(9) 



-V3a 2 b* (3|a| 2 - 2) a (3|a| 2 — 1) & ^3a*b 2 
V3ab* 2 -(3|a| 2 -l)6* (3|a| 2 -2)a* V3a* 2 b 
-b* 3 V3a*b* 2 -V3a* 2 b* a* 3 

We present only a quarter of the transition matrix for S — 2. The rest can be found by 
using above described symmetry properties. 



(8) 



/ a 4 V4a 3 b 



(4|a| 2 -3)a 2 y/6(2\a 



U 2 



V6a 2 b 2 Vlab 3 b A \ 

2 -l)ab (4|a| 2 -l)6 2 
6|a| 4 -6|a| 2 + 1 



(10) 



III. DENSITY MATRIX AND BLOCH TENSORS. 

When the random magnetic field acts onto the spin, the system must be described by the 
density matrix p. By definition it is a (2S + 1) x (2S + 1) Hermitian matrix with the trace 
equal to 1. It satisfies the standard equation of motion: 



Any Hermitian matrix with the trace equal to 1 can be represented as a sum: 

P = 2sTT J + s • S+bik (s i s k + S k Si - lS tk S{S + 1)) + 
■■■jisy9hi2--i2s {ShSi2---Si2s + a ^ permutations — all traces) 



(11) 



(12) 



5 



If the Hamiltonian is the generator of the rotation each term in equation (|12p corresponds 
to an irreducible representation and evolves independently. We will call symmetric tensors 
9iki 9iM ■■■9iii 2 ...i 2 s the Bloch tensors in analogy with the Bloch vector g well known from 
Bloch theory of the nuclear spin motion. Any trace of such a tensor must be equal to zero. 
The Hamiltonian (J2J) generates following equations of motion for the Bloch tensors: 

g = — b x g; g ik = —EilmblQmk ~ £klmbl9im', ■■■ (13) 

All these equations are independent and have obvious integrals of motion: 



g 2 = const; g\ k = const; g\ kl = const; ... (14) 

Thus, the density matrix of a spin S in an external time-dependent magnetic field has 
25 conserving values. It is convenient to represent the Bloch tensors by their complex 
component with definite projection to the z— axis. We will denote such components of 
a tensor of the rank s as g s ,m- The corresponding tensor operators composed from the 
symmetrized products of 2s components of the spin S operators are denoted T s 5 m . They 
can be constructed from the senior operator of this representation T^ s = 2~ S / 2 S^ with 
S± = S x ± iS y by recurrent commutations with the operator S-. 

rpS _ \ 

J(s + m + l)(s — m) 

The operators T s 5 m are polynomials of the standard spin operators S± and S z . They are 
operator analogs of spherical harmonics. Several lowest such operators are presented in 
Appendix. We show below relations between the Cartesian components of the tensor gi 1 ,„i s 
and its components g S)Tn for several values of s: 



' s,m+l 



(15) 



91, ±i = ~^(9 X± Wv) ' 9i,o = 9z (16) 

92, ±2 = -J= (gxx ~ 9yy ± 2tg xv ) ] 0j,±l = ~J= (9xz ± %9yz) ] 92,0 = 9zz (17) 
^3,±3 = -J= {9 xxx i 3ig X xy — %9xyy T Wyyy) '1 ■■■ (18) 



The general rule for writing the s, ±m— component via its Cartesian counterparts is the 
same as for the product (x ± iy) m z s ~ m . The Hamiltonian (j2J) in terms of components 



with definite projections reads (note that b± = b x ± ib y ; S± = S x ± iS y ): 

H = -b z S z - l -(b + S- + bS + ) (19) 
Equations ()13|) in terms of the components with definite z— projections read: 

9s,m = -imb z g Sjm + ~yJ (s + m)(s - m + l)b + g s , m „i + ^(s -m)(s + m+ l)b^g S)in+1 (20) 
and the conservation laws are: 

s 

\9s, m \ 2 = const (21) 

m=—s 

IV. FAST NOISE IN TWO-LEVEL SYSTEM 

In this section we consider only spin 1/2 or, equivalently a two- level system. We assume 
that the magnetic field can be separated into regular and random parts: 

b(t)=b r (t)+Tj(t) (22) 
where b r (t) = zb z t + xb x and r](t) is the Gaussian noise determined by its correlators: 

(Vi{t)vk(t')) = fik{t-t') (23) 

We assume that the correlators (|2l?|) decay after a characteristic time difference r n and 
that this correlation time is much less than the characteristic time of the LZ process tlz- 
However, the noise must be slow enough to avoid the direct transitions between the levels 
when the interlevel distance approaches its saturation or characteristic value u far from the 
crossing point. Thus, the noise correlation time r n must satisfy a following inequalities: 

LO- 1 < T n < T LZ (24) 

The spectral width of noise is l/r n . The noise produces transitions during the interval of 
time t acc = l/(b z r n ), after which the current LZ frequency becomes larger than the noise 
spectral width. We will call this interval the accumulation time and assume that it is much 
larger than other characteristic time intervals r n and tlz- 
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We first solve an auxiliary problem in which b x = and transitions are mediated by 
noise only. Such a problem for a special shape of correlators (f xx = J 2 exp (— the 
remaining components of the correlation tensor are zero) was solved earlier by Kayanuma^ 
and studied numerically by Nishino et alM. In our work 16 we have generalized and simplified 
the Kayanuma solution. Here we reproduce our solution^ 6 - in a modified form convenient 
for extension to higher spins. We also obtain new results calculating the fluctuation of the 
density matrix, or equivalently the Bloch vector. Equations for the component of the Bloch 
vector in this case are: 

g z = (i/ V2) (r] + g_ - r/_g + ) ; g± = =p (b z t + rj^ g± + (if V2)r]±g z (25) 

Solving equation for g±, we find: 



/ ib z t 2 r \ 
9±(t) = g±(-oo)exp =F— ^— =F « J r) z (t')dt' \ + 



i/V2) 



exp 



ri±(t')g z (t')dt' 



(26) 
(27) 



Let us first consider the case of complete initial decoherence: g±(— oo) = 0. Then, plugging 
equation (j2Tj) into the first equation (j2SJ), we find a separate equation for g z : 



g z = -(1/2) / exp 



ib z (t 2 -t' 2 ) 



- i / Vz(t")dt h 



r)+(t)r)-(t')g z (t')dt' + c.c. 



(28) 



Let us average equation (J28j) over the ensemble of the random noise. For such averaging it is 
important that the noise correlation time r n is much shorter than the time t acc necessary for 
a substantial variation of (g z ). This fact allows to represent the average (r] + (t)ri-(t')g z (t > )) 
approximately as a product: (r) + (t)ri_(t')g z (t')) = (r] + (t)r]_(t')) (g z (t')). More accurately one 
should incorporate the fluctuations of g z . In the leading approximation they are determined 
by the same equation (}28|) as follows: 

ib z (t\-t\ 



8g z = -(1/2) J dtt J exp 



-oo — oo 



(v + (ti)v-(h) - (v+(ti)v-(h))) (g z (h)) dt 2 + c.c 

(29) 

We ignore rj z (this approximation will be justified by the next step). Let substitute this 
additional term into equation (j2~Hj) and perform averaging over the gaussian random field 
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T]. According to the Wick's rule, it is reduced to all possible pairings. In our case the only 
possible pairing is (r)+(t)r)-(t2)) (V-ityv+iti))- Such a pairing limits the integration by the 
interval t — r n < t 2 < t\ < t' < t. Thus the contribution of the fluctuational term differs 
by an additional factor ~ T n /t acc << 1 from the principal contribution from (g z ). These 
arguments represent a shortened version of the original arguments by Kayanuma^ and are 
akin to the Abrikosov-Gor'kov theory of static disordered alloys 20 . 

Using the fact that the decay of the correlator (rj + (t)r]-(t')) limits effectively the inte- 
gration over time by an interval t — r n < t' < t, we can prove that the contribution of the 
noise component r] z in the exponent in equation (|28j) can be neglected. To estimate this 
contribution we assume that r\ z is statistically independent from other components. Then 



f Vz (t")dt" 



. The 
evel of 



the averaging over r] z results in the Debye- Waller factor exp 

argument of this exponent can be estimated as (r/ 2 ) r 2 . It is small provided the 
noise (rjl) is much smaller than r~ 2 . The noise correlators are even functions of the time 
difference. Therefore, expanding linearly the time argument of the exponent in the same 



equation 



b z (t 2 -t' 2 ) 



b z t(t — t'), one can transform the integral-differential equation (J28J) into 



an ordinary differential equation for (g z ): 

(g z ) = —F (b z t) (g z ) 
where F(Q) is the Fourier-transform of the function F(t) = fxxij) + fyy{ T ) 

oo 

F(n) = J F(r)cosfirdr 

— oo 

Equation (jBUjl has a simple solution: 



(30) 



(31) 



(9z(t)) = (&(-oo))exp 
At t — > +oo the asymptotic value of (g z ) is: 

(g z (+oo)) = (g z (-oo))exp(- 



F (b z f) dt' 



(32) 
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F(0) 



(33) 



Note that what matters for the LZ transition is the average quadratic fluctuation of non- 
diagonal noise at any moment F(0) = (r/l + Vy") i n contrast to a standard characteristic of 
the white noise which would be F(0). Indeed, commonly white noise correlator is introduced 
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as (rj(t)r](t')) = j5(t — t'). The only characteristic of the noise in this approach is 7 = 
J!f£o (vityvft)) dt' . An interesting feature of the asymptotic formula (jHHj) is its independence 
on the noise correlation time r n . However, it should be kept in mind that this asymptotic is 
valid only at time t t acc = (b z r ri ^ 

Returning to the solution (|27|) for g±, we see that the term 
g±( — 00) exp ^ = F^~ "F * / Vz(t')dtj > omitted at substitution in the first equation 
f!25jl and the ensemble averaging, vanishes if the z— component of the noise is statistically 
independent from others. Now let us perform a similar procedure solving first equation for 
g z and then substituting the solution into equations for g±. In the same approximation 
equation for averages of these components of the Bloch vector reads: 



{g ± ) = -\P (b z t) (g ± ) 



(34) 



The first term in equation (|27j) after the averaging turns into zero at any finite t. Indeed 

/ / 00 \ 2 \ 

the Debye- Waller factor which appears in this case is exp 



asymptotics of (g±) at t — > +00 is: 

(9±(+oo)) = exp 



f Vz (t")dt" 



(g±(-°°)) 



vrF(O) 



0. The 



(35) 



Note that the symbols g±(±oo) denote the coefficients at exp 

(t 1 ^)- Thus, the noise 

asymptotically tends to reduce the average components of the Bloch vector, i.e. to establish 
equipopulation of the levels and to destroy the coherence. However, during the time interval 
of the order of t acc the average Bloch vector can oscillate. 

For the considered problem it is possible to find exactly the fluctuations of the Bloch 
vector. Indeed, they are given by a standard formula: 



[(5 g y) = {^) - ( g y (se) 

Due to the symmetry of the problem, g 2 is a conserving value. Therefore, its average 
coincides with itself and is determined by initial conditions. The average value of the vector 
(g) was calculated above. Thus, we find an asymptotic value of the fluctuations: 



t=+oo 



^(-00) (1 - exp (-29)) + g 2 x (-oo) + g 2 (-oo) (1 - exp (-9)) (37) 
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The values of average square fluctuation can be also written for any moment of time. The 
fluctuations are strong, i.e. their magnitudes are of the same order as the average values of 
the Bloch vector components unless 9 is very small. An important property of the noise is 
that in fluctuations it mixes diagonal and non-diagonal elements of the density matrix, i.e 
all three components of the Bloch vector. 

Let us proceed to the solution of a more general problem retaining non-zero x— component 
of the regular field b x . Such a solution becomes possible due to separation of times 16 : we 
can neglect the effect of b x beyond the LZ time interval t^z near the crossing point. On the 
other hand, we can neglect the effect of the noise inside this and even much larger interval of 
time, since its characteristic time is t acc ^> tlz- Thus, the problem is separated into 3 parts: 
in two intervals (—00, —to) and (t , +00) we can use the solution of our auxiliary problem, 
i.e. to take in account only transitions caused by the noise; in the interval (— to, to) the 
Landau-Zener solution is valid. We need only to match them properly. The time interval to 
must satisfy a strong inequality: tlz *C to "C t acc . As we have shown before, the average 
components g z and g± evolve separately under the action of the noise. Therefore, at the 
moment of time —to these components are: 

(9-A-to)) = exp(-0/2)<7 z (-oo); (g±(-t )) = exp(-0/4)</±(-oo) (38) 

The action of the LZ transition matrix (£Q) with the matrix elements (J7J) can be transferred 
onto the vector g by using spin-1 matrix (JBJ). Thus, the average components of the Bloch 
vector at t — to are as follows, 

(g z (t )) = V2ab* (g+(-to)) + (2|a| 2 - 1) (g z (-t )) + V2a*&<0_(-t„)) (39) 
(<7+(*o)> = ^ (<7+Ho)> - V2ab (g z (-t )} - b 2 (g-(-t )) (40) 
(S_(t )) = -b* 2 (g+(-to)) - V2a*b* (g z (-t )) + a* 2 (<7_Ho)> (41) 

where a and b are given by equations (J7J). The transition from +t to +00 is provided by 
the same diagonal transition matrix (|38jh which was already used for the transition from 
—00 to — to- 

(<fc(+oo)) = exp(-0/2)(fc(+t o ); (<7±(+oo)) = exp(-0/4)<7 ± (+* o ) (42) 
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Collecting together equations (j38H42|) . we find the final result: 



(g z (+oo)} = v^exp (-36/4) (ab*g+(-oo) + a*&<?_(-oo)) + exp (-6) (2|a| 2 - l)<? z (-c^3) 
(g + {+oo)) = exp (-6/2) (aV(-oo) - &V(-°°)) - exp (-30/4) V2abg z (-oo) (44) 
(g-(+oo)) = exp (-6/2) (-b* 2 g + (—oo) + a* 2 g_(-oo)j - exp (—30/4) \/2a*&*<k(— oo)(45) 

Let us analyze first the case of complete decoherence at t = — oo, i.e. g±(— oo) = 0. Then 
equations (|43H45|) look much simpler: 



(g z (+oo)) = exp(-0)(2|a| 2 -l)^(-oo) (46) 
(g+(+oo)) = -ex V (-36/A)V2abg z (-^) (47) 

Equation (|4T)|) shows that, in the absence of initial coherence the population difference can 
only decrease after the transition. The noise only strengthens this tendency. However, 
the initial coherence if exists can increase the difference of population to the value larger 
than initial. Equation ()47|) demonstrates an inverse process: the coherence can appear after 
the transition even if it was absent in the initial state. It is straightforward to derive the 
transition probability from equation (|46jl: 



1 



1 



P^ 2 = - [1 - exp(-6)(2\a\ 2 - l)j = - [1 + exp(-0) - 2exp(-0 - 2irY)\ (48) 

This formula was obtained in our previous work 16 . At zero noise intensity (6 = 0) this 
result turns into the Landau-Zener transition probability. At very big noise (6 = oo) the 
probability is equal to 1/2. The probability is the same for transitions 1 — > 2 and 2 — > 1. 
This symmetry does not look strange keeping in mind that we considered the classical noise, 
which produces the induced transitions. It vanishes as soon as the quantum nature of the 
noise and spontaneous transitions are taken into account. 

For completeness we present a formula for quadratic fluctuations at t = oo in the case of 
complete initial decoherence: 

(g 2 (+oo)) - (g(+oo)) 2 = |^(-oo)| 2 [l - e- 2d - 4 (e- 2 ^ 2 - e^) (e^ 2 - e~ 2d )] (49) 

The fluctuations vanish at 6 = and reach their maximum value equal to \g z (—oo)\ 2 at 
6 = oo. 
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V. FAST NOISE AT A MULTILEVEL CROSSING. 



We consider only the case of a Zeeman multiplet placed into a varying magnetic field. It 
is described by equations (|2(J|) . The averaging of them over the fast noise is performed by 
two steps as it was done in the previous section. First we neglect the transitions produced 
by the regular part of the magnetic field and take in account only the transitions produced 
by the random field. This approach is correct outside the time interval r^z near the avoided 
crossing point t = 0. Assuming for simplification complete initial decoherence, we find for 
the average diagonal matrix elements following equations: 



9s,o( +oc )) = ex P 



s (s + l) e 



, M-oo)) (50) 

Note that these relations do not contain S explicitly, the time evolution depends on s only. 
In the course of derivation of equation (jBTIj) we ignored not only the fluctuations of gf Q , but 
also all higher components except of gg±i- Their contributions to the main component 
have the same order of magnitude ~ T n /t acc as the contribution of fluctuations. 

The generalization of equation (|5U|) to the higher projections m, i.e. to the coherence 
factors requires some care. The truncated system of equations (|2()jl for zero initial values of 
all components of the s-tensor except of gf m , reads 



i>s,m 
Js,m—1 
ys,m+l 



-imb z tg s s m + i (X s ,mV+9s, m -i + X s ,-mV~9 
-i(m-l) b z tg\ 3 x + i\ m r]_g s 



s 

s.m+l I ' 



-i (m + 1) b z tgl 



,m+l 



+ iK 



iV+9s 



(51) 
(52) 
(53) 



where A s>m = y (s — m)(s + m + 1). To deal with slow- varying average values the fast 
oscillating exponent should be eliminated. In order to do that we introduce slow variables 
9s,m 9s,m 6Xp (^fm) . Further we consider only slow amplitudes g s , m and omit the sign 
tilde. After elimination of the values gs m ±i and averaging, we find a following equation for 
q s : 



q S 



(s{s + 1) - m 2 ) F(b z t) + mG(b z t) 



(54) 
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where F(Q) is defined by equation (J31)) and G(Q) is defined as a sine Fourier-transform: 



(^x(t)^(O) - Vy{i~)Vx{fy) smttrdr 



(55) 



Thus, the time dependence of the average (gf m ) is defined as: 



9s,m(t) 



exp 



s(s + 1) - m 2 ) F(b z t') + mG(b z t')} dt' (<?f m (-oo 



(56) 



Its asymptotic value at t — > +oo does not contain the sine Fourier-transform: 



exp 



(s(s + I) - m 2 ) 6 



-oo. 



(57) 



At s = l,m = 1 this result coincides with 

Now we proceed to our main problem including both regular and random transverse 
magnetic fields. We will employ the separation of their action in time proven in the previous 
section. To avoid lengthy formulas our consideration will be restricted to the case of complete 
initial decoherence. Then the only non-zero components of the Bloch tensors are gf . Their 
evolution is described by three independent factors, two of them originating from the noise 
and the central factor being the generalized Landau- Zener-Hioe matrix element: 

s(s+ 1) 



exp 



-0 



P^(2e 



(58) 



where P s °'°(x) is the Jacobi polynomial. The average values of the Bloch tensors components 
with m/0 vanish as a result of averaging over the random phases in the initial state. To 
find the transition probabilities Pj->j> it is necessary to put all the diagonal elements of the 
density matrix except of pjj equal to zero in the initial state: 



1 



2S 



1 s=l 



-oo 



jki 



(59) 



and find from these equations the initial values (7* (— oo). Then the transition probabilities 
are: 



2S 



- 1 s=l 



-OO 



(60) 



2S + 1 ' ^" s ' uv '"' \- a *>Jjj'~ r \ 2 
where g^ (+oo) are defined by equation (J58)l. In the following tables we demonstrate the 
results of calculations according to this algorithm for the values of S 1, 3/2 and 2 (the results 
for S = 1/2 are given in the previous section by equation. (jjHj) )- 
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+1 



+1 \ + ~Ei + \E 2 



3 



1 ~ i E 2 



I + 1^2 



-1 



(61) 



3/2 



1/2 



3/2 


i 

■i 


+ 


9Ei 
20 


+ 


£ 2 
4 


+ 


#3 
20 












1/2 


i 

■i 


+ 


20 




E-> 
4 




20 


1 + 


20 


+ 


E% 
20 


i 9E 3 
~ 20 


-1/2 


i 

i 




20 




4 


+ 


3£ 3 
20 


1 

4 


Ei 
20 


+ 


#2 

20 


9E 3 
20 


-3/2 


i 

4 




9-Ei 
20 


+ 


K 2 
4 




20 













-1/2 -3/2 



(62) 



A J 








2 


2 


i 

5 


+ 


2Ei 
5 


, 2£ 2 _i_ E a 
~ 7 ~ 10 


1 


1 

5 


+ 


Bi 
5 


£2 £3 
7 5 





1 

5 




2£ 2 
7 


1 3E 4 
35 


-1 


1 

5 




El 
5 


£2 1 E 3 
7 " r 5 


-2 


1 

5 


+ 


2Si 
5 


, 2£ 2 1 £3 

7 10 



+ 



£4 

70 



2g 4 1 
35 5 



Mi 
10 



1 1 E2 
5 7 



2E 4 1 
35 5 

, Ei 
~ 70 



10 



Ez 4. 

14 " r 

12E 4 
35 

la _ 

14 



2£ 3 
5 



2^3 

5 



8iJ 4 
35 



8E4 
35 



-1 -2 



1 1 2^2 1 18E 4 
5 7 35 



(63) 



100 



Here we denoted E s = exp (-^^-&), i.e. E 1 = e - e ; E 2 = e" 3e ; E 3 = e" 6e ; E 4 = e 
Unfilled sites in the table can be easily restored using the time reversal symmetry: = 

Pj'-^j = P-j-^-j'- 

In Figl we compare our formulas with the numerical solution of the Schodinger equation 
for the spin 5=1, which convincingly confirms our analytical results. 

The quadratic fluctuations of the values g~ m for fixed S and s are calculated as it was 
done in the previous section: 



* E k 



s 



E k 



s 



t=— 00 



In the case of complete initial decoherence employing equation ([580. we find: 



(64) 



* E k 



s 



9s A 



-00 



t=+oo 



I -El (P?>° (2e- 2 ^ 2 -1 



(65) 
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FIG. 1: The final probability to find a spin-1 system in the state with S z = — 1, or 1 when the 
initial state is S z = 1 as function of the noise amplitude. The Hamiltonian is H = tS z + rj x S x 
where (rix(ti)r] x (t2)) = J 2 e _A ' il_ * 2 '. Discrete points correspond to results of numerical simulations 
with averaging over 200 different noise realizations and A = 125. Lines correspond to analitical 
predictions of (|61|). 

VI. NOISE AT ADIABATIC LEVEL CROSSING. 

In this section we consider a situation of adiabatically changing levels, i.e. b z -C but 
still the noise is supposed to be sufficiently fast r n <C {b z )~ l l 2 . The adiabatic Landau-Zener 
transitions in coupled two level systems have been proposed as a candidate for implementa- 
tion of quantum gates in quantum computing 21 . It is important to understand the influence 
of noise on such a q-bit manipulations in order to control the error propagation in quantum 
gate circuit. We do not specify the relationship between b x and r n . This study is motivated 
by experiments with nanomagnets^*^ in which inequality \fh z T n <C 1 can be easily realized 
since the sweeping rate of the applied magnetic field can be made arbitrarily small. However, 



16 



the noise may be not fast enough to compete with the tunnelling amplitude b x . The nu- 
clear bath correlation time is in the range of r n ~ 1ms, whereas the measured values of the 
tunnelling amplitude for known nanomagnets range between 10~ 10 — 10~ 3 K, or equivalently 
10 — 10 8 Hz. In a part of this interval b x r n has a rather large value. 

We consider in this section only 2-level systems. A natural approach to this problem is the 

transfer to the adiabatic set of states, i.e. to the eigenstates of the instantaneous regular part 

I a x {t) 

of the Hamiltonian (fI9|l. Let denote this time-dependent eigenvectors as a(t) = 



and b(t) 



a 2 (t) 



-a 2 (t) 
ai(t) 



, where ai (t) = J eJ ^; a 2 (t) = J £ -^ and e(t) = Jb 2 z + ¥ x , 



b z t). The unitary matrix of transformation to the adiabatic set reads: 



U(t) 



a\{t)I — ia 2 {t)a y (66) 



a>i(t) a 2 (t) 
y -a 2 (t) ai(t) 

where a a (a = x, y, z) are Pauli matrices. In the new basis the total Hamiltonian acquires 
the following form: 

H = e{t)a z + UrjaU- 1 = e{t)a z + rj'a (67) 
The transformation of the random field components is: 

/ / 2 2\ i o b z b x / t b x b z (ro\ 

Vx = i a i - a 2>Vx + 2aia 2 ri z = —i] x + —i] z ; i] y = rj y ; r\ z = -~Vx + —Vz (68) 

In this form the Hamiltonian (JBTj) essentially coincides with the Hamiltonian of the auxiliary 
problem (zero regular transverse field) for the two- level system (see Section ITVj) . The essen- 
tial difference is first, that the effective regular external field is not linear in time; instead 
it is equal to e{t)] second, the correlators of effective noise rj' now depend not only on the 
time difference, but also on time itself due to the time-dependent transformation (jSHJ): 



<»/«(*K(0> 



bitt'( Vx (t) Vx (t')) + b 2 x M*MO> + 



e(t)e{t' 

b z b x t (rj x {t)rj z (t')) + b z b x t' {r}z{t)rix{t')) 



(69) 
(70) 



Still the noise correlation time is small in comparison to the characteristic time of variaton 
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for the adiabatic energy e{t). Employing the same approximation as in Section HVl we arrive 
at a similar equation of motion for average in the case of complete initial decoherence: 



(9'M) = ~F\t) (g' z (t)) 



(71) 



where 



F\t) = fyy (6(f)) + 



+ b z b x t 



Kt'f xx (e(t)) + Kf zz (e(t)) 



(72) 
(73) 



In the last equation the hats symbolize Fourier-transforms of corresponding correlators. As 
before, we can find the average value (g z (t)) at arbitrary moment of time. Asymptotically 
at t — > +oo we find: 



(g , z(+ 00 )) = ex P 



b 2 z t 2 L (e(t)) + e\t)f yy ( e (t)) + b\f zz ( e (t)) 



(9'A-oo)) (74) 



The characteristic time after which the correlators in equation (J71j) become very small and 
decay rapidly is determined by approximate equation e(t)r n « 1. If r n <^ b^ 1 , then this 
characteristic time coincides with the accumulation time t acc = (b z T n ^j defined in Section 
IIV| terms proportional to b x are negligibly small and we return to the result (|4*3*|) with 
\a\ = 1, or equivalently to (jHSJ) of the Section ITVl In the opposite case r n 3> b~ x the value 
e(t) exceeds r^ 1 at any moment of time t. Therefore, all correlators are small and the value 
of exponent in (|74|) is close to 1. It means that practically no transition proceeds due to 
the noise between adiabatic states. Thus, equation (f?4*j) carries most interesting information 
when r n ~ b^ 1 . An interesting feature of the transition probability is that the z-component 
of noise can produce transition between adiabatic states. This happens because the latter 
rotate with time. Note that ^-component of noise is irrelevant if r n <C b~ x . 



VII. CONCLUSIONS 



Motivated by synthesis and magnetic measurements of cubic nanomagnets, we developed 
a theory which allows to find the transition probabilities between the states of the Zeeman 
multiplet in the presence of the regular time-dependent and random magnetic field (noise). 
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The solution of this problem occurs to be possible since the evolution matrix for the quantum 
problem is a rotation matrix acting in a spin S representation. The density matrix can 
be expanded into a linear superposition of irreducible tensor operators. The coefficients 
at this operators related to one of the irreducible representations (Bloch tensors) evolve 
independently on others. Thus, the initial problem in the space of dimensionality (2S + 1) x 
(2S + 1) is reduced to 2S separate problems in the linear spaces of dimensionality from 1 to 
2S. 

The second key observation is that, for the fast noise, the transitions due to the noise and 
those due to the regular part of magnetic field are separated in time. This fact allows to solve 
the problems for regular field and for the noise in the absence of the regular non-diagonal field 
separately and then match them. An interesting conclusion of our theory is that, in contrast 
to usual statistical calculations with the white noise, in which only the Fourier-component of 
the noise correlation function with zero frequency matters, the transition probabilities in the 
Landau-Zener problem depend only on the average square of the random field amplitude. 
We were able to find asymptotically exact analytical results for the probabilities. From them 
we concluded, that, in the absence of initial coherence, the average values of the diagonal 
components of the Bloch tensors (with zero projection onto the direction of sweeping field) 
decrease monotonously with time due to the noise. It means that the population differences 
in average can only decrease after the transitions. However, they can grow if there is a 
coherence in the initial state and non-diagonal components of the Bloch tensors are not 
zero. Due to high symmetry the considered system has 2S additional integrals of motion: 
traces of the square of each Bloch tensor. Thus, the increase of the population differences 
proceeds at the expense of the non-diagonal components, i.e. coherence amplitudes and vice 
versa.. 

The same conservation laws enabled us to find exactly the fluctuations of the Bloch's 
tensors, in particular the fluctuations of the transition probabilities in the genuine Landau- 
Zener problem. They are of the same order of magnitude that the average values. 

The noise in our theory is considered as the classical random field. To incorporate the 
quantum properties of noise is an interesting and challenging problem. The second unsolved 
problem is to study the correlation of the Bloch tensors at different moments of time. 
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IX. APPENDIX 

Here we present several lowest operator spherical harmonics Tg . In contrast to scalar 
spherical harmonics they depend on 3 parameters, two of them are integers (s and m), 
whereas S accepts integer and half-integer values. The simplest nontrivial harmonics are 
vectors with s = 1. They are: 



rpS C . rpS C 

Next we demonstrate second order tensorial harmonics: 



(75) 



Tr 



2,±2 



1 Q2 . rpS 



2 •-=■=' -2,±1 — 
The third rank harmonics read: 



- (S±S Z + S Z S±) ; T 2fi 



si--s(s + i) 



(76) 



rpS 

J 3,±3 



2372^±' ^3,±2 — {s%Sz + S±S Z S± + S Z S±^J 



T i,±i - \l 12 



S Z S+ + S Z S + S Z + s+s. 



3S(S+1) - 1 



3,0 



3 3g(g + l)-l 



(77) 
(78) 
(79) 



For the fourth rank harmonics we find: 



J 4,±4 — 



V 



■4,±3 



rj^ ^S±S Z + S±S Z S± + S ± S Z S± + S Z S: 



2 5/2 



(80) 
(81) 
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rpS 

A -t-9 — 


V7r 
12 L 


rpS 

1 4,±1 — 


y/7 

2 5/2 


rpS 

J 4,0 ~~ 


4 



S Z S± + S Z S±S Z S± + S Z S±S Z + S±S^S± + S±S Z S±S Z + S±Sz — kS^. (82) 
+ S 2 S ± ,S 2 + S 2 S ± ^ 2 + 5±5^ - k (5 2 5± + S±5 2 )] (83) 
(5* _ K S 2 Z + A) (84) 



where we have introduced notations: 



6S(S+ 1) -5 _ 3g(g + l)[g(5+l)-2] 
K " 7 ; " 35^ 



(85) 
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